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Colloidal solutions posses a wide range of time and length scales, so that it is unfeasible to keep 
track of all of them within a single simulation. As a consequence some form of coarse-graining 
must be applied. In this work we use the Multi-Particle Collision Dynamics scheme. We describe a 
particular implementation of no-slip boundary conditions upon a solid surface, capable of providing 
correct forces on the solid bypassing the calculation of the velocity profile or the stress tensor in 
the fluid near the surface. As an application we measure the friction on a spherical particle, when 
it is placed in a bulk fluid and when it is confined in a slit. We show that the implementation of 
the no-slip boundary conditions leads to an enhanced Enskog friction, which can be understood 
analytically. Because of the long-range nature of hydrodynamic interactions, the Stokes friction 
obtained from the simulations is sensitive of the simulation box size. We address this topic for the 
slit geometry, showing that that the dependence on the system size differs very much from what is 
expected in a 3D system, where periodic boundary conditions are used in all directions. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

The existence of a huge range of time and length-scales, 
spanning between mesoscopic colloidal particles and mi- 
croscopic solvent particles, constitutes a severe problem 
in numerical simulations. Hybrid schemes have been de- 
veloped, in which different coarse-grained approaches are 
used to describe solvent molecules and colloids. Promi- 
nent in the class of solvent models is Multi-Particle Col- 
lision Dynamics (MPCD), originally proposed by Male- 
vanets and Kapral jlj, which has proved to be very ef- 
fective in simulating Newtonian fluids out of equilibrium. 
For a review see 0]. 

The interaction between solvent molecules and colloids 
may be described in several ways, mimicking the various 
boundary conditions (BC) in use in continuum descrip- 
tions. The molecular origins of these boundary condi- 
tions can be very complex [3, |1] , but for most colloidal 
applications, simple no-slip BC are sufficient for model- 
ing experimental conditions. In this paper we discuss a 
way to impose such boundary conditions in MPCD sim- 
ulations of liquids containing dissolved colloids. 

MPCD simulations consist of two alternating steps: 
the streaming step (where Newton's equations of motions 
for non-interacting particles are solved) and the collision 
step (where the fluid is coarse-grained). The no-slip BC 
involve both steps. In particular, during the collision 
step, virtual particles (VP) are inserted in the regions of 
the MPCD box occupied by walls or colloids. This idea 
has already been used in p-Q. In this paper, we will 
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provide additional insight in this method, analyzing not 
only the velocity profiles nearby solid planar walls, but 
also discussing the contributions of the VP to the forces 
exerted by the solvent on the solid boundaries. Such 
forces determine the drag forces experienced by the col- 
loids and thereby regulate the colloids' dynamics. A re- 
cent interesting work discussing these techniques in the 
implementation of mixed stick-slip boundary conditions 
can be found in [§]. 

The second topic we treat is the role of finite-size ef- 
fects (FSE) of the simulation box. In molecular dynamics 
simulations it is common practice to use periodic BC to 
mimic an infinite system. However, due to the long-range 
nature of the hydrodynamic interactions, quantities such 
as the Stokes friction can still be affected by the system 
size even when equilibrium properties are well reproduced 
[l0| . The use of periodic BC implies a set of images of 
the colloidal particle, so that they all together form a 
periodic grid. The problem of slow flow through a 3D 
periodic array of spheres has been treated by , while 
in [l^l FSE are studied for a chain of beads in a solvent. 
In both works, periodic BC are assumed along all three 
directions. 

When confining walls are present in one direction, and 
periodic BC in the others (such as in the simulation of 
a slit geometry), the problem of the FSE has been ad- 
dressed only recently in jl2i | for slip BC, using molecular 
dynamics simulations, and in [l^ for no-slip BC, using 
continuum models. Our results, obtained with coarse- 
grained simulations, are qualitatively in agreement with 
such previous works for the friction in the direction paral- 
lel to the walls. Furthermore we also discuss the friction 
in the direction perpendicular to the walls, which is miss- 
ing in jl^ . [13I ] . For both the parallel and perpendicular 
friction strong differences emerge with respect to the case 
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of a cubic box with periodic BC in all directions. 

The structure of this paper is as follows: in section 
the simulation technique is introduced, with particular 
emphasis on the use of the VP during the collision step; 
in section IIIIl we validate the model through the study 
of the Poiseuille flow in a slit, for which the theory is 
known. In particular we calculate the forces exerted by 
the solvent on the solid walls through two independent 
methods, once by using the VP and once using the stress 
tensor for a MPCD fluid near planar boundaries. In sec- 
tion IIV[ we compute the friction on a single sphere in 
bulk. Because of the VP, the local Brownian friction is 
increased, and we provide an Enskog-like model to pre- 
dict such an effect. In section |V] the colloidal particle is 
confined in a slit, and we measure the friction as a func- 
tion of the lateral size of the walls, while keeping fixed 
the separation between the walls. Final remarks and ob- 
servations are in section IVII 

II. SIMULATION TECHNIQUE 

In the present application, Multi-Particle Collision Dy- 
namics (MPCD) is a hybrid simulation scheme, in which 
a coarse-grained approach is used to describe the solvent 
variables, while an atomistic description is adopted for 
the solvent-solute and for the solute-solute interactions. 
The dynamics of the system is made up of two steps: 
streaming and collision. In the streaming step, the posi- 
tion and velocity of each particle is propagated for a time 
St by solving Newton's equations of motion. In the col- 
lision step the fluid is subdivided into cubic cells of side 
a. Then a stochastic rotation of the particles velocities, 
relative to the center of mass motion of the relevant cell, 
is performed according to the formula: 

V^{t)=U + n{v,{t)-u}, (1) 

where u is the mean velocity of the particles within a cell 
and O is a matrix which rotates velocities by a fixed an- 
gle (in this work a = tt/2) around a randomly oriented 
axis. Through the stochastic rotation of the velocities, 
the solvent particles can efficiently exchange momentum 
without introducing direct forces between them during 
the streaming step. As the collision step conserves mass, 
linear momentum and energy, the correct hydrod yna mic 
behavior is obtained on the mesoscopic scale [H Tl^. as 
long as a shifted-grid procedure is included in order to en- 
force Galilean invariance [l^. When colloids are present, 
Newton's equations of motion are solved also for them 
. Special care must be used if no-slip boundary con- 
ditions (BC) are applied on the colloidal surface and on 
the confining walls. In this case, in fact, we must couple 
colloids and walls to the solvent during the collision step 
too. This is achieved by means of virtual particles (VP). 
The implementation of the no-slip BC is described in the 
following. 

Streaming step: when a MPCD particle crosses the col- 
loid (or wall) surface, it is moved back to the impact 



position. Then a new velocity is extracted from the fol- 
lowing distributions for the tangent (vt) and the normal 
components (v„) of the velocity, with respect to the sur- 
face velocity: 

pM = T-rpexp {-mvl/2kBT) , (2) 



p{vt) — y^m/2TrkBTexp (— muj /2fcsr) , (3) 

where m is the mass of the solvent particle, fc^ Boltz- 
mann's constant, and T the temperature of the system. 
Once the velocity has been updated, the particle is dis- 
placed for the remaining part of the integration time step. 
Collision step: VP are inserted randomly in those parts 
of the system which are physically occupied by the colloid 
or by the walls (in sufficiently thick layers behind the in- 
terfacial positions). The VP density matches the MPCD 
solvent density 7, while their velocities are obtained 
from a Maxwell-Boltzmann distribution, whose average 
velocity is equal to the velocity of the colloid surface or 
to the velocity of the walls, and the temperature is the 
same as in the solvent. According to their coordinates, 
VP are sorted into the grid cells. During the collision 
step, the average velocity of the center of mass of the cell 
is computed as 

riMPCD rivp 

nMPCD + nyp 

where umpcd and nyp respectively are the number of 
MPCD particles and the number of VP belonging to the 
same cell. Finally velocities of both MPCD and VP be- 
longing to the same cell are rotated according to the rule 
given in Eq.((T]). 

Due to the exchange of momentum between the sol- 
vent and the colloidal particle, the force exerted upon 
the latter may be expressed as: F = fs + fc, where fs 
and fc are the forces during the streaming step and the 
collision step, respectively. The former can be calculated 
as: 

f^ = -^E^P'''''^''' (5) 

i—l 

where Q is the number of MPCD particles which have 
crossed the surface of the colloid between two collision 
steps and /\pMPCD jg -j-j^g change of momentum of the 
«-particle of the solvent, which has been scattered by the 
colloid. The force exerted during the collision step is: 

= ^E^Pr, (6) 

where q is the total number of virtual particles which 
belong to a tagged colloid; and Apy^ is the change 
of momentum of the i- virtual particle during the col- 
lision step. If walls with no-slip BC are present, we can 
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FIG. 1: (Colors online) Velocity profile across a slit of width L^. 
Simulations are performed in a box of sides Lx = Ly = = 64a. 
The external applied field is = 10~^/oz. The inset shows a 

magnification of the profile near the wall at a; = 0. 
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FIG. 2: (colors online) Average force exerted by the fluid upon each 
wall of a slit of width Lx = 64a, versus time. The total expected 
value is f<^^v = 16/o; fs and /c are the force exerted during the 
streaming and the collision steps. 



also measure the force exerted by the solvent upon them 
through Eqs. In such a case Q is the number of 

solvent particles that have crossed the surface of a wall 
and q is the number of VP belonging to the same wall. 



III. POISEUILLE FLOW THROUGH A SLIT 

In this section we study the flow of the MPCD solvent 
through a slit under the influence of an external uniform 
force F'^^*, which is oriented parallel to the walls. Walls 
are placed at a; = and x = and the lateral sides 
of the walls are Ly = L^- In our simulations we choose 
the solvent density equal to 7 = 5, and the interval be- 
tween two collision steps equal to St — O.Uq. The time 
is in units of Iq = a{m/kBTy^'^ , where m is the mass 
of the solvent particle, ks the Boltzmann constant and 
T the temperature. Hereafter we assume that fcsT = 1 
and m — 1. A stationary parabolic velocity profile is 
expected to form for an incompressible fluid. Simulation 
results for the case F^^* = with /o — ksT/a 

are plotted in Fig. ([T|). When the no-slip BC are imple- 
mented only in the streaming step, using the stochastic 
reflections according to Eqs. ([21131), the slope of the ve- 
locity profile (solid line) just near the wall is different 
from the expected one (dotted line); moreover large slip- 
page near the wall persists. When VP are included into 
the collision step, the Vz{x) profile (dashed line) repro- 
duces quite well the Navier-Stokes solution. The use of 
VP clearly ameliorates the simulation results with no-slip 
BC as expected 0-13 • We next ask what are the conse- 
quences of the various changes for the force exerted by 
the fluid upon the wall? From a balance of external body 
and wall forces in a stationary state, the expected value 
of the force on each wall is /'^^p = -fF"''* kfLyL^ [13, 
which in our case is f'^^'^ = I6/0. 



In Fig. (21) we see that the total force converges to the ex- 
pected value of /'^^P = I6/0, at which point a stationary 
state is reached. In the stationary state the contribution 
to the total force due to the collision step (/c) is rela- 
tively large. This happens especially when the mean free 
path of the MPCD fluid is small, in which case momen- 
tum is transported mainly via collisions rather than via 
diffusion. 

For the slit geometry it is easy to adapt the general 
expression for the stress tensor of the MPCD fluid with 
periodic BC [l^ [l^ in order to measure the force just 
near the wall. The force exerted along the z direction 
because of collisional exchange of momentum along the 
X direction is: 

= iE^eA-... (7) 

In this expression we have used the same notation as 
adopted in [15], where a full description of the MPCD 
stress tensor is provided. The meaning of the nota- 
tion is as follows: A^|'^ = ^j,^{t + St) - ^|'^(i + St), 
Avj^z = Vj,zit + St) — Vj,z{t + St); Vj^z and Vj,z are the z 
component of the velocities after and before the collision 
step of particle j. ^j{t) = {^j,x,^j,y,£.j,z) are the coordi- 
nates of the cell of the unshifted grid containing the par- 
ticle j at time t. Similarly are the coordinates of the 
cell of the shifted grid containing the j particle at time t. 
The sum in Eq. ([7]) applies only to the p particles of the 
MPCD solvent belonging to the cells which overlap one 
of the walls. As shown in Fig. ([2]), the expressions of the 
forces according to Eq. © and ([7]) are in very good agree- 
ment. However, if one needs the force exerted by the fluid 
upon a spherical colloid or an irregularly shaped object, 
it is important to take into account the local curvature 
of the colloidal surface. In such a case the generalization 
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FIG. 3: (Colors online) Running integral of the autocorrelation 
functions of the constraint force. The off-diagonal terms of Eq. (ISj, 
which are not shown here, go quickly to zero. 



of the stress tensor for the MPCD fluid is not easy, while 
the implementation of Eq. ^ is straightforward. 



IV. FRICTION ON A SPHERE IN THE BULK 

We now turn to the translational friction on a sphere 
which is not free to rotate. In this respect we consider 
the case of a single sphere of radius Rcoi = 4^cl, which is 
kept at a fixed position by applying a constraint force 
This condition is similar to the case of very mas- 
sive particles embedded in solvent as studied in [9|. We 
use a cubic simulation box whose side is L = 32a. We 
are primarily interested in the autocorrelation function 
of J-'^ because this can be connected, via a Green-Kubo 
relation, to the translational friction tensor : 



lim - — 

t->oo ks 



T 



dT{T^^{to + T)J^I{to)). (8) 



where a, /3 G {x, y, z}. 

Because of symmetry, the friction does not depend 
on the cartesian direction along which we measure, see 
Fig. ([3]). If we consider a particular direction a, the long 
time limit of provides the total friction ^ upon the 
sphere, for which there are essentially two contributions: 
one coming from the local Brownian collisions with the 
particles of the fluid {£,e), while the other is due to the 
long-range hydrodynamic interactions (^s). A simple 
empirical formula says that the hydrodynamic and the 
Brownian friction should be added in parallel [l^, [2l| in 
order to obtain the total friction: 



1/e = i/^E + i/^s 



(9) 



Both ^ and may be calculated using our simulations. 
^ can be read from the long-time limit of the integral in 



Eq. (8), while corresponds to the height of the short- 
time peak [22| . The hydrodynamic term can be ex- 
tracted by inverting Eq. ^ . In the present case we obtain 

= 1088 ± 4 and .^J"" = 270 ± 30. 
The hydrodynamic friction can be theoretically estimated 
from the drag force Fd on a colloid with no-slip BC in 
an infinite fluid medium, according to the Stokes law 
Fd — —£,'^Voo = —&TTrjRcoiVoo, where v^o is the flow field 
at large distances. As the Stokes friction depends on the 
long-range hydrodynamic effects, it can be substantially 
affected by the finite-size L of the simulation box. De- 
tailed calculations in [10] suggest that to lowest order in 
Rcoi/L, the correction is given by = £.^l3 (Rcoi/L), 
with 



p{Rcoi/L) « (1 - 2.837 (i?eoi/i))" 



(10) 



Taking into account the finite-size effects, the value for 
the Stokes friction is = 291, and the simulations result 
is in good agreement with this. 

We now turn to the Brownian friction S^e , due to local 
collisions with the solvent. A suitable starting point to 
calculate this term is the Enskog-Boltzmann theory for 
a dense gas. This model takes into account only two- 
bodies collisions and successive collisions are supposed 
to be uncorrelated. In other words no influence of the 
local disturbance, induced by the tagged colloidal particle 
on its surroundings, is included. The expression for the 
translational friction in a bath of particles of mass m is 
the following [H: 



8 ( 2nkBTmM \^'^ 2 l + 2x 
" 3 m + M ) ^^'°'TT^ 



(11) 



where x = ^ I^^R^coi ~ 2/5 is the gyration ratio for a 
sphere. According to this model S^e ^ 700, while from 
simulations we obtain a much larger value. The discrep- 
ancy is a consequence of the changes made to the collision 
step, as the exchange of momentum between the MPCD 
particles and the colloid, via VP, changes the effective 
number of local Brownian collisions. 

In order to evaluate the fluctuations in the constraint 
force due to the VP, let us focus on one collision 
cell which partly overlaps with a colloid and contains p 
MPCD particles and q VP, with velocities and w^, 
respectively. Let us suppose that the colloid moves with 
velocity V through the ideal gas bath, then we know that 
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V 



m 

2\ 



3fcsT 



(12) 
(13) 



The collision step itself may be written as (primes indi- 
cate velocities after the collision): 



u = 




i=l 

w\ = V, + (O - I) (v, - u) 
W- = Vifj + (Jl - I) (Wj - u) 



(14) 

(15) 
(16) 
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where ft represents a rotation of an angle a around a ran- 
dom axis, I is the unit matrix, and u the average velocity 
of the cell center of mass. The total change of momen- 
tum of the VP, i.e. their contribution to the change of 
momentum of the colloid (or wall) is equal to: 



AP 



VP 



= m(ri-I) 



p 



p + q 



q 



(17) 



We define the force exerted on the colloid by the cell 
under consideration as F^^ = AP^^St and consider its 
contribution to the Enskog friction matrix: 



H(i) = (3 / (F(r)F(0))dr 



- B {F{0)F{0)) 6t 



3=1 



(F(j(50F(0)) St. 



(18) 



Only the first term will be non-zero because velocities 
at different times are uncorrelated. When calculating 
this term, we assume < ViVj >— ^^^ISij, < WiWj >= 
^^ISij and < ViWj >— 0, which turns H = H(0) into a 
diagonal matrix. Averaging over all possible orientations 
of the rotation axis we obtain for the diagonal elements: 



e 2 m pq 

ivp = ^(1 - cosa) — — — . 
6 ot p + q 



(19) 



Since velocities in different cells are uncorrelated we may 
simply add the contributions of all cells overlapping with 
the colloid, obtaining: 



(tot 



-(1 — cosa)-^ > . 

3' 'Stc Pk + qk 

cells A; 



(20) 



Performing the sum over cells analytically is rather 
complicated, if not impossible, so we have decided to 
evaluate Eq. pO)) numerically, during the simulation it- 
self. Including the correction so obtained, the predicted 
short time friction becomes: = + ^vp — 1084, 
which is in very good agreement with the simulation 
valued;"" = 1088 ± 4. 



V. FRICTION ON A SPHERE INSIDE A SLIT 

Recent developments in the field of micro- and nanoflu- 
idics have led to a renewed interest in molecular hydro- 
dynamics phenomena near a solid surface. It has been 
shown theoretically (23 - [27l | and experimentally (28l - [3lj | 
that the hydrodynamic interactions with fiat walls can 



slow down the motion of colloidal particles substantially, 
and that such effect also depends on the direction of the 
motion of the particle (for example parallel or perpen- 
dicular to the walls). We therefore apply our method to 
study the friction on a sphere inside a slit. In this case, as 
in the bulk system, we do not take into account the role 
played by the angular momentum, so that the present 
discussion is relevant only for purely translating spheres 
in a slit. The walls are set at z = and z — Lz while we 
use periodic BC along the x and y directions. The length 
of the walls' sides is Lx = Ly = L. We have chosen this 
particular configuration because we have access to ana- 
lytical solutions of the Navier-Stokes equations for the 
parallel motion and numerical solutions for the per- 
pendicular motion as well (27l | . Such solutions have been 
provided under the assumption that the lateral size (L) 
of the walls is infinitely large. 

Our aim is to study how the friction depends on the 
lateral width L of the cell when the walls separation is 
fixed. Besides establishing finite size effects for the usual 
simulation boxes, this topic is also relevant for the study 
of experimental periodic arrangements of particles, such 
as trains and grids of particles under confinement. In 
both cases a periodic array of colloidal particles exists, 
exerting strong hydrodynamic interactions on each other. 
Simulations have been performed keeping the colloidal 
particle of Sec lIVI in a fixed position in the mid-plane 
of the slit. From the autocorrelation of the constraint 
force we obtain the friction coefficients (as explained in 
section lIVp in the directions parallel and perpendicular 
to the walls. Under the assumption of small Reynolds 
numbers, it is possible to express the effects of the walls 
by means of a correction factor r to the Stokes friction 
in an unbounded system: 



it 



(21) 



Similarly we describe finite-size effects by a second cor- 
rection factor B: 



/-_L P-L„ L /-oo 

^2 = P ^ ?S ■ 



In the following we will obtain as 



/-^ isim 1 ^1 



(22) 



(23) 



Similar expressions hold for the parallel frictions. 

Analytical expressions for the parallel and perpendic- 
ular frictions in a slit geometry are not available, but 
analytical expression for a single wall do exist [1^. The 
correction for the parallel friction due to the presence of 
a single wall is: 

rll = [1 - 1.004(i?eo;A) + 0.418(i?co//z)'+ 

+0.21(i?,o,/z)4 - 0.169(i?eoi/z)'] , (24) 
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FIG. 4: Finite-size effects corrections for the friction upon a spher- 
ical particle placed in the mid-plane of a slit geometry. The lateral 
size of the square walls is L, while the wall separation is Lz = 32a. 
Circles: simulation results. Solid line: the corrective terms for a 
cubic box with periodic BC in all three directions (Eg llOI I. 



FIG. 5: (Colors online) Finite-size effects for L > 20 in comparison 
with two different fit functions. Red lines are based on a function 
of the type as in Eq. mUl l: black lines are based upon the function 
in Eq. II27I I. Dotted lines are for the perpendicular friction, dashed 
lines for the parallel one. 



whereas the correction for the perpendicular friction is: 

n{n + 1 



(2n- l)(2n + 3) 



= - sinh(a) 

2 sinh(2n + 1)q; -|- (2n -|- 1) sinh 2q; 
4 sinh^ (n -f 0.5)a - (2n 1)2 sinh^ (a) 
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is very close to its images, a wake effect is established 
which reduces the friction along the lines connecting the 
colloid and its nearest images. A similar trend can be 
deduced from [H, [l^ . The situation depicted in these 
works is slightly different from ours. They for example 
study the parallel friction as a function of once Ly 
and Lz are fixed, while in the present simulations we al- 
(25) ways use a square periodic space {L^ ~ Ly). However if. 



where a = cosh^^ (z/Rcoi) and z is the distance of the 
particle from the wall. As a first approximation, we use 
linear superposition theory, according to which we can 
add the contributions coming from each wall as if they 
behave independently from one another: 



' w2 



1, 



(26) 



the subscripts wl and w2 refer to the two walls. We have 
subtracted 1, in order to avoid counting the bulk contri- 
bution to the friction twice. This is a mere consequence 
of the definition adopted for the corrective terms. 

In case the walls are very far apart, this approximation 
works fairly well, otherwise it tends to overestimate the 
combined effects of the two walls. In particular, for the 
simulation parameters we have used here {Rcoi — 4a and 
Lz = 32a), the friction obtained with the superposition 
model appears to be overestimated by about 15% with 
respect to the numerical solution of the Navier-Stokes 
equations in the slit as given in [27]. In the rest of this 
section we will use the results provided by Eq. (j26|) but 
corrected according to ^3 ■ 

The values of /3^'ll are plotted in Fig. (|4]) as a function 
of lateral size L. 

For values of L < 24a, the behavior of /j" and 
is very different: with decreasing L, Z?" decreases while 
increases. For the parallel motion we expect that, 
when the lateral size L is very small and the particle 



for each profile shown in [12l [13[ , we pick up the points 
corresponding to L^ = Ly, it is seen that the friction de- 
creases as L diminishes. On the contrary when L > 24a, 
the parallel and perpendicular frictions behave similarly: 
they decay slowly to the case of a slit whose walls are in- 
finitely large. It is possible that the screening point where 
we can observe different behaviors between the parallel 
and perpendicular friction is for L ~ L^. 

When L = Lz, the simulation box is cubic, so we can 
compare our results with those provided by the correc- 
tion factors in Eq. pUj) plotted in Fig. as a solid line. 
We see that there exists a significant difference which is 
due just to the presence of the confining walls: the lat- 
ter screen even the effects due to the use of periodic BC. 
However, once the wall separation Lz is fixed, the cor- 
rection factor decays slower when compared to Eq. (fTU]) 
where periodic BC are used in all directions (Fig. 

For ease of use, we have tried to fit the data for L > 20 
with a function similar to that in Eq. (jlOp (see red lines in 
Fig. (O): the behavior seems good for the perpendicular 
friction but it deteriorates for the parallel friction. We 
have also tried a fit with the function 



f{L) = l + ylln(l + i?,„j/L), 



(27) 



in which A is a parameter whose value depends on the 
direction we consider: for the parallel friction a" ~ 2.08 
while for the perpendicular motion ^ 3.02. Eq. ((27|) 
seems to work better than Eq. PH)) . For example the 
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parallel friction is better captured (see black lines in 
Fig. (O). We do not have, however, indications whether 
the parameter A is a function of the wall separation 
and of the particle position z. To clarify this point fur- 
ther simulations are necessary, to be compared with the 
solutions of the Navier-Stokes equations for the slit for 
different particle positions. 

In summary, our results suggest that: 

• There are two different regimes depending on 
whether the lateral size of the walls is smaller or 
larger than the wall separation. 

• For a cubic box, the finite size effects in the slit 
geometry are less than in a system with periodic 
BC in all three directions. 

• Once Lz is fixed and L > L^, the correction fac- 
tors decay slowly towards unity, possibly as a 
logarithmic function of L. 

VI. CONCLUSIONS 

In conclusion we have described a protocol for the 
MPCD simulation technique, which provides a satisfac- 
tory treatment of the no-slip boundary conditions (BC) 
and the correct evaluation of the force exerted by the sol- 
vent on solid surfaces. Such a protocol is based upon the 
use of virtual particles during the collision-step, a method 
already used in literature and which appears suited to 
study wall-liquid interfaces and very massive particles 



embedded in the solvent, or particles kept fixed by an 
external constraint force as in the present work. We pay 
particular attention to the calculation of the friction ten- 
sor for a non rotating sphere. 

We have shown how the no-slip BC modifies the local 
Brownian friction on a spherical colloid, and how such ef- 
fects can be evaluated through an Enskog-like treatment 
of the VP. We have also studied the friction on a parti- 
cle embedded in a slit, analyzing the dependence of the 
simulation results on the use of periodic boundary condi- 
tions along the sides of the walls. When the lateral size 
of the walls is very small (less than the wall separation), 
strong coupling effects between a particle and its images 
are observed. Moreover, the parallel and the perpendic- 
ular friction show an antithetic behavior: as L decreases, 
the parallel friction decreases too, while the perpendicu- 
lar friction increases. When L > L^, the corrective terms 
for the parallel and perpendicular friction are much more 
similar to each other and they both decrease to unity as L 
increases, approaching more and more the case of an ideal 
slit made of two infinitely large walls. When L = Lz, we 
find that the finite-size effects (FSE) in the slit geometry 
are less than for a cubic box with periodic boundary con- 
ditions in all directions. Moreover, once Lz is fixed and 
L > Lz, the FSE in the slit appear to be slowly varying 
with the system size L, possibly according to a logarith- 
mic function of Rcoi/L. 

Further studies will concern whether and how the FSE 
in a slit depend on the walls' separation Lz and on the 
particle distance z. 
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